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Abstract 

We present the evaluation of the two-loop 0(g 4 M t 2 /M^) effects in the partial widths of 
the Z° boson in the MS scheme and in two different implementations of the on-shell 
scheme. We observe a clear reduction of the scheme dependence of the predictions. The 
renormalization procedure and the Heavy Top Expansion employed in the 0(g 4 Mf /M^) 
calculations are illustrated in some detail and intermediate results are provided. We 
discuss the implication of our results on the constraints for the Higgs mass making use 
of simple interpolating formulas. We find that precision data give M H < 285 GeV at 
95% CL. taking into account the theory uncertainty. Including also the information from 
direct search experiments we obtain a 95% upper bound M H < 345 GeV. 



1 Introduction 



The impressive amount of data collected at LEP, SLC, and the Tevatron allows for a 
very stringent test of the Standard Model (SM), at the level of the per mille accuracy. 
This experimental situation has required a significant effort from the theory side in order 
to provide predictions for the electroweak observables with comparable precision. We 
recall that such a level of accuracy tests the quantum structure of the theory, thereby 
providing information about physics at mass scales that are not directly accessible at our 
experimental facilities. 

After the discovery of the top quark at Tevatron the only mass scale still unknown in 
the SM is the one related to the Higgs boson. As the sensitivity of the various precision 
observables to the Higgs mass, M H , is approximately logarithmic, the extraction of the 
relevant information is more difficult and delicate than in the case of the top and requires a 
careful consideration of higher order effects. In the recent past it was recognized that 
in order to obtain accurate constraints on M H from global fits to the experimental data, 
the theoretical predictions of the observables most sensitive to M H should also include 
two-loop top effects enhanced by factors (M 2 /M^,) n (n — 1, 2). This observation called for 
the study of the 0(g 4 M 2 /M^) corrections to the M w — M z interdependence, presented 
in Ref. |J, and to the effective mixing parameter sin 2 6 l $ H- Moreover, in Ref. § the 
analysis of the 0(g 4 M 2 /M^) contributions to M w and sin 2 d^fl was carried out in the 
MS formulation of Ref. f| as well as in two different implementations of the on-shell (OS) 
scheme ||. It was shown that the inclusion of 0(g 4 M 2 /M^,) effects sharply reduces the 
scheme and residual scale dependence of the theoretical predictions for M w and sin 2 O^f}, 
and suggests a comparable reduction of the overall theoretical uncertainty |7|[|. 

The aim of this paper is to go further in the study of the 0(g 4 M 2 /M^) corrections. 
Our purpose is actually manifold: 

• First, we present the result of the calculation of the partial widths of the Z°, Tf, 
for / ^ b, at 0(g 4 M?/Ml) in both the MS and OS schemes. As the case of 
the b quark involves a completely new class of vertex diagrams and would require a 
significant dedicated effort, it will not be considered here. We also study the scheme 
dependence of the results for the leptonic width, 1^, finding it significantly reduced 
by the inclusion of the new contributions. 

• We take the opportunity of the new result on the partial widths to describe some 
important aspects common to all 0(g A M? /M^) calculations which were not dis- 
cussed in the preceding publications, in particular concerning the renormalization 
procedure and the Heavy Top Expansion (HTE) technique used in the calculations. 
We also provide some intermediate results and derivations in appendix A. 

• Finally, we discuss some of the consequences of our results, in view of the present 
experimental situation. Following Ref. J7|, we provide elementary interpolating 
formulas in different schemes for T^. For completeness, we also report and update 
the analogous formulae for M w and sin 2 O^fl. 
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A few topics relative to the 0(g A M^ /M^) calculations were already covered in some 
detail in preceding related publications, and we will not dwell on them further. For 
instance, the impact of the new calculations on the indirect determination of M H is quite 
important, and this has been extensively discussed in Refs. Although we do not 

repeat this discussion here, we will use the interpolating formulae for M w , sin 2 6^*, and 
T^, together with the latest experimental data, to obtain two new estimates for M H . One is 
based on information coming only from precision measurements, while the other includes 
also the results of direct search experiments, using the method developed in Ref. ||10|| . 

The paper is organized as follows. In the next section we discuss some issues concerning 
the renormalization procedure we have followed. Section || deals with the technical matter 
of the HTE. In section |] we present the analysis of the corrections to the Z partial 
widths to 0(g 4 M 2 /M^). Section 5 contains the interpolating formulae and our estimates 
for M H . Finally, we present a brief summary. We also include two appendices. In 
appendix A we provide the two-loop W and Z° self-energies to 0(g 4 M 2 /M^), which may 
be useful in different contexts. Appendix B contains instead the explicit expressions of 
the 0(g 4 M 2 /M^) corrections to the various Z partial width form factors. 

2 Renormalization 

In this section we discuss some aspects of the 0(g 4 M 2 /M^,) renormalization for processes 
involving four light fermions as external states. The main difference between the one and 
two-loop cases is that in the latter we also have to consider the renormalization of the 
unphysical sector. In the following, we first discuss the renormalization of the unphysical 
sector in a general way. We then explain why the use of OS masses for the vector bosons, 
rather than MS mass parameters, is particularly useful in our calculations and illustrate 
different options for the renormalization of their unphysical counterparts. 

As is well known, there are constraints that link the longitudinal component of the 
vector bosons, the corresponding pseudo-Goldstone bosons, and the Faddeev-Popov (FP) 
ghosts. In higher order calculations, it is convenient to choose a renormalization procedure 
that automatically respects the Slavnov- Taylor Identities (STI) which are induced by the 
local gauge invariance of the original Lagrangian before spontaneous symmetry breaking. 
According to the organization of the calculation, it is possible to use different procedures 
that respect the STI's and are particularly convenient in order, for instance, to minimize 
the number of diagrams to be considered. Of course, physical amplitudes are independent 
of the chosen procedure, and this can be used as an additional check of the calculation. 
For the problem at hand, the discussion can be kept at the one-loop level, and we can 
limit ourselves to contributions proportional to powers of M t . In the following, we discuss 
only the W sector; the case of the Z° sector can be treated in a completely analogous 
manner, because the 7 — Z° mixing does not contribute at the order of our calculation. 

If we split the unrenormalized W polarization tensor into transverse and longitudinal 
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parts 

n^(g) = (gr - ^r) A ww (g 2 ) + n^.(g 2 ), (1) 

and we denote by U. W( f, the proper two-point function for the mixing between the W and 
its pseudo-Goldstone boson, 0, and by the self-energy of the latter, we obtain the 
following STI's in a general gauge (see for ex. UTTjpl 

e fe ^(g 2 ) + M z U w<j) (q 2 )} + G(q 2 ) = 



Q 2 



Kw(r) + 2M w U w<p (q 2 ) + M 2 W U^(q 2 ) + M 2 W T = Q, (2) 



where T represents the tadpole contribution and G(q 2 ) is a term involving charged FP 
ghosts as external fields. Among other things, G(q 2 ) contains the two-point function of 
the charged ghosts; as the ghost do not couple directly to fermions, at the one-loop level 
G(q 2 ) is completely independent of the top quark. 

We see that the first two terms of the second STI in Eqs. (0) vanish at q 2 = 0, and 
that therefore 11^,(0) is equal to the tadpole contribution. This uncovers the connection 
between the renormalization of the Goldstone boson mass and the one of the tadpole. In 
our calculation we employ the usual tadpole renormalization that minimizes the effective 
potential and consists in removing all tadpole graphs. This choice implies the subtraction 



of II^(O) from the two-point function of the pseudo-Goldstone boson Hl3|.|Tll . 

Concerning the longitudinal component of the two-point function of the W, we find it 
convenient to renormalize it in the same way as the transverse part, the one which is re- 
lated to the mass counterterm. Therefore, the choice of the renormalized mass parameters 
employed in our analysis is going to affect the structure of the finite part of the coun- 
terterm in the unphysical sector. In practice, the choice is between the MS and OS mass 
definitions. For the other two-point functions appearing in Eq. (0) different options are 
possible, which all respect the STI's, and are equivalent at the level of physical amplitudes. 
They correspond to different ways of renormalizing the gauge-fixing parameters. 

A first and very natural option (heavy mass procedure) consists in subtracting from 
each term in Eqs. (0) the contributions proportional to powers of M t . In case an MS 
prescription for the vector boson masses is adopted, the subtraction applies only to the 
divergent part. This procedure is equivalent to a Taylor expansion in the external mo- 
mentum and obviously respects the identities, as the unrenormalized self-energies do. We 
recall that I1^ W (0) = A ww (0) ~ Re A WW [M 2 7 ) at the order we are interested in, so that 
the OS mass counterterm corresponds to the first term of the Taylor expansion. There is 
a counterterm for the W — <fi transition as well as for the <f> two-point function, but the 
FP ghost self-energy needs not be renormalized, as it is independent of the top quark. 
In order to satisfy the second of Eqs. (Q) at any q 2 , a wave function renormalization 
for the <f> field is needed. On the other hand, <f> appear only inside the loops in our 
calculations and the use of is indeed not necessary. Therefore, we simply demand that 
the second of Eqs. @ be satisfied at q 2 = M^. 

A second possibility consists in assigning no counterterm to the W — <p transition; this 
decreases the number of counterterms and simplifies their implementation. In fact, in the 
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Figure 1: Two-loop vertex diagrams containing the top. Wavy (dashed) lines represent vector 
(scalar) bosons. 

't Hooft-Feynman gauge it corresponds to renormalizing the masses of the vector boson 
and of the associated scalar boson and ghost in the same way (bare gauge fixing), apart 
from the supplementary subtraction at q 2 = 0, corresponding to the tadpole contribution, 
that we have discussed above. It therefore amounts to fix = 5M^ + T and 5M 2 = 
5M^ (c^ 1 are the FP ghosts) and is the closest to the naive parameter renormalization. 
Also this choice verifies both Eqs. @ at q 2 = M^, leaving room for a further arbitrary 
wave-function renormalization, which again can be avoided altogether. 

We have explicitly verified that these two renormalization options (heavy mass and 
bare gauge fixing) are equivalent at the level of physical amplitudes in the cases at hand. 
Finally, we should mention that another possibility is the one presented in [15|, which 
corresponds to the use of renormalized gauge-fixing. 

For what concerns the three-point functions, we just recall that the renormalization of 
those couplings that contain masses is fixed by the Ward identities that link the Yukawa 
and bosonic coupling counterterm to the gauge coupling and mass renormalization. 

From a technical point of view, the choice of OS masses as renormalized parameters, 
in particular for the W and Z° vector bosons, has some important advantages in the 
calculations considered here. The processes we are interested in always involve light 
external fermions and are mediated at lowest level (to excellent approximation) by gauge 
bosons only. In these cases, the two-loop box and vertex diagrams containing the top 
quark consist of a top-bottom insertion on a bosonic line belonging to a one-loop vertex 
(Fig.|I]a,b) or box[|, plus vertices that involve a W — transition, like the one shown in 
Fig.|l]c. Superficially, all these diagrams are either convergent or logarithmically divergent. 
However, they contain a quadratic subdivergence associated with the top-bottom loop 
insertions. If we employ the OS mass definition for the mass of the vector bosons, all 
the boxes and vertices involving insertions on vector boson lines can be neglected in the 
calculation because the associated mass counterterm graphs remove completely their M 2 
dependence. Of course, this would not happen in the case of a pure MS subtraction, as 
the remaining finite part would still be proportional to M 2 . 

The use of OS vector boson masses therefore absorbs a whole class of 0(g 4 M 2 / 'M^) 
diagrams and we have decided to present always our results in terms of OS masses for 
all the physical particles, also when we work in the MS scheme. Indeed, this choice 
has some clear advantages even in the case the gauge couplings are renormalized in the 
MS scheme |5],[nj: (i) it employs as mass parameters the physically measured masses 

X A top loop vertex insertion on a trilinear coupling does not provide M t 2 terms by dimensional reasons. 
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without introducing large radiative corrections in the coupling counterterms; (ii) it absorbs 
the large radiative corrections that appear in the one-loop relation between OS and MS 
masses. Incidentally, one should also note that the conventional OS mass definition, 
identified by the counterterm Re A WW (M^), does not present any gauge ambiguity at the 
order of our calculations; as stressed in ||17|| , this is not the case at 0(g A ) and beyond. 

The heavy mass procedure that assigns a counterterm to the W — <ft transition has 
an additional advantage: it allows us to neglect systematically all two-loop vertices and 
boxes in our calculations, including the diagrams of Fig.lc. In this case the only two- 
loop diagrams that we need to consider are two-point functions of the vector bosons: 
this is the approach followed in the present paper. If instead we work with bare gauge 
fixing, it should be apparent that only diagrams of the kind in Fig.^(a,b) will be removed 
at 0(g A Mf /M^) after renormalization. The diagram of Fig.|](c) has to be explicitly 
calculated. However, as could be expected, its evaluation at 0(g 4 M t 2 /M^) involves only 
self-energy integrals and therefore no additional difficulty is introduced. This approach 
was followed in Refs. 0,|3| and the detailed evaluation of the diagram of Fig.|I](c) was 
discussed in Ref. [p~8|1f|. 



3 Heavy top expansion of two-loop diagrams 

In the preceding section we have seen that in the calculations at hand the only two-loop 
diagrams that need to be considered are two-point functions of the vector bosons, including 
the 7 — Z° mixing for the calculation of the effective sine and the photon self-energy at 
vanishing momentum for the renormalization of the electric charge ||. In this section we 
explain how we have evaluated these self-energy diagrams. We start motivating the choice 
of the next-to-leading order (NLO) HTE instead of a numerical computation, proceed to 
consider the classes of diagrams that need to be computed, and conclude discussing the 
quality of the approximation. 



There are arguments |22|] that two-loop self-energy integrals with arbitrary masses 
and momentum transfer cannot be expressed in terms of known elementary functions like 
polylogarithms. Although an analytic solution seems in the general case impossible, there 



exist several methods that permit an efficient numerical evaluation [p3| . On the other 
hand, compact analytic results are available in a number of special cases, which apply 
to a variety of physical situations like small and large momentum expansions. The clear 
advantage of the analytic evaluation in a situation where at least one of the parameters, 
namely M H , is largely unconstrained, is that the results can be readily implemented in 



independently of the renormalization of the unphysical sector, vertex diagrams of the kind shown in 
Fig.[l](c) can be excluded by the choice of a non-linear gauge fixing; for example the r yW ± (/> ' 1 coupling, 
which is characteristic of the gauge, can be avoided, so that the photon separately couples to W and 
charged Goldstone bosons jl9) . Consequent modifications [^0) in the Faddeev-Popov sector of the SM do 
not affect the calculation, as ghosts do not couple to fcrmions. In the case of the photon vertex, also the 
background field gauge can be used to the same effect. Alternatively, it is possible to fix the gauge in 
such a way that the Z°W ± (j) T vertex is avoided pi|. 



5 




(a) 00 



,o. 



/A/VWWW 



.O. 



(e) 

Figure 2: Two-loop topologies for vector boson self-energies. Dashed (solid) lines represent 
virtual bosons (fermions). 

fitting routines. Furthermore, if one chooses to perform a numerical calculation at the 
two-loop level in the electroweak theory without any approximation, several challenging 
problems arise, like the very large number of diagrams, the implementation of the complete 
two- loop renormalization of the SM, etc. Therefore, following the conclusions of []T],0], our 
strategy has been to use a heavy mass expansion up to 0(g A M 2 / 'M^). 

We have performed the calculations in the 't Hooft-Feynman gauge, where = M c = 
M w . This choice proves to be convenient also in the evaluation of the integrals, as it limits 
the number of different mass scales occurring in each diagram. Concerning the bottom 
mass, we have set it to zero from the start. We have neglected all flavor violation. All 
two-loop self-energy diagrams containing the top have to be considered. The topologies 
are displayed in Fig.^J The subset of diagrams containing only scalar bosons in addition to 



the top and the bottom quarks is responsible for the leading 0(g 4 M A /M^) correction (24 

We distinguish between three classes of diagrams: (i) the ones for which a naive 
Taylor expansion in q 2 can be used because no light mass threshold is involved; (ii) those 
containing at least one bottom quark threshold, which we regard as zero mass threshold; 
(hi) those with at least one physical threshold of the order of the W mass (e.g. it can 
be 4M^, M^, [M w + M z ) 2 ). Clearly, the two last classes of diagrams overlap (see for 
ex. the diagram in Fig. 0(c) where the fermion loop consists of two bottom and one top 
lines and the bosonic lines represent two W bosons). With respect to this classification 
of diagrams, the fact that the Higgs boson mass is a priori arbitrary implies that we have 
to consider two scenarios: (a) the case in which M H is light compared to the top mass, 
namely M H « M w and therefore M H <C M t \ (b) the case in which M H is considered 
heavy compared to M w , and the heavy mass expansion is implemented leaving the ratio 
between M H and M t finite and arbitrary. For most diagrams containing the Higgs, the 
HTE provides different expressions in the two cases. This explains why in app. A the 
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expressions for the W and Z° self-energies are given in the two regimes. 

In the case of class (i) diagrams we adopt a Taylor expansion first in the small external 
momentum and then in the light masses keeping only terms proportional to powers of M t . 
The Taylor expansion of the two-loop integrals is viable for values of q 2 up to the first 
physical threshold. Since in our calculations q 2 = z or 0, it follows that q 2 <C M 2 
and this procedure is justified for class (i). Typical diagrams belonging to this class are 




Figure 3: Two-loop master self-energy diagram. 

those of Fig|J(a) where the fermion in the loop is a top quark. Referring to the two-loop 
master diagram of Fig.[| any diagram containing a heavy particle (top and possibly Higgs 
boson) in place 1 and 2 can be safely included in this class. In practice, the q 2 expansion 
is implemented using (k and m generically indicate one momentum of integration and one 
internal mass) 



1 



EC- 1 ) 



(2k-q + q 



2\n 



(3) 



(k + q) 2 — m 2 ^ (k 2 — m?) n+l 

under the n-dimensional integral. For arbitrary internal masses, the two-loop integrals at 
zero momentum transfer can be expressed as a combination of logarithms and dilogarithms 
25] . Therefore, we reduce the integrand to scalar vacuum integrals and for the evaluation 



of the integrals with higher powers of the propagators we follow the iterative method of 



first paper in [25 



For the diagrams of class (ii) the Taylor expansion is bound to fail because it cannot 
describe the physical threshold at q 2 = 0. The typical diagram has bottom at place 1 
and 4 and top at place 2 and 3 in Fig. [| The simple replacement of Eq. (^) when m = 
leads to increasingly severe infrared divergences for increasing n. Heavy mass or small 
momentum expansions in these cases are therefore asymptotic expansions, about which a 
rich technical literature exists [E(|. An elegant algorithm to implement the asymptotic 



expansion in dimensional regularization has been developed in Ref. [2?J and applied to 
this class of diagrams in P7fl . The technique we have used, however, is based on the 
repeated application of the identity 



k\ + 2fci 



(fci + k 2 ) 2 - M 2 k 2 - M 2 [(kt + k 2 ) 2 - M 2 ] {k 2 - M 2 ) 



(4) 



where M is a heavy mass and are the integration momenta. The first term on the 
r.h.s. of Eq. (Q) leads to a disconnected integral (product of two one-loop integrals) that 
can be evaluated exactly, i.e. for arbitrary q 2 , and contains lng 2 terms that describe the 
physical threshold. The second term, instead, leads to a two-loop integral with improved 
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infrared convergence in the ki integration and improved ultraviolet convergence in the h\ 
integration. Therefore, if, for example, the original integral is IR logarithmically divergent 
in the ki integration when q —>■ 0, the corresponding two-loop integral associated with 
the second term in Eq. (Q) evaluated at q 2 = gives a finite result that differs from 
the complete result for arbitrary q 2 by terms 0(q 2 /M 2 \n(q 2 /M 2 )). In case of integrals 
more seriously IR divergent, an iterative use of Eq. gives an approximate evaluation, 
up to terms suppressed by the heavy mass, through a combination of integrals (a set of 
disconnected ones plus two-loop integrals evaluated at q 2 = 0) whose analytic expressions 
are known. Thus, one can combine the expansion (|3|) with the repeated use of Eq. (^) 
to obtain the asymptotic expansion of the original integral at the desired order. Indeed, 
at each step of the expansion the more serious IR divergent behavior introduced by (§) 
can be balanced by the convergent IR factor from the second term of Eq. @, while the 
part associated with the first term in Eq. (Q) can be always evaluated exactly. An explicit 
example of the application of this method has been provided in ||18||. We have checked 



that our results coincide with those of Ref. |27] in all the cases relevant to the present 
calculation. 

The class (iii) diagrams differs from those of class (ii) because the threshold is not 
at q 2 = and can in fact be well above the q 2 at which the integrals are evaluated. 
Nevertheless, the Taylor expansion cannot be employed here because we are interested in 
the HTE and not in the small momentum expansion. Therefore, in order to extract all 
the terms enhanced by powers of M t , one has to use the same procedure adopted for case 
(ii), reducing the original integral to a sum of disconnected ones, to be evaluated exactly 
keeping the light masses (M w , M z and possibly M H ) and the q 2 , plus two- loop connected 
integrals with q 2 = which have to expanded in powers of M 2 . The disconnected diagrams 
are responsible for the terms containing the functions Bo [q 2 , m\ , m 2 ,] in the two-loop self- 
energies of app. A. Also in this case we have checked our results against the relevant 



existing literature [28 



A crucial question is now whether the approximation based on the NLO HTE described 
above is reliable. The second term 0(g 4 M 2 / M^) of the two- loop HTE seems to be quite 
important wrt the first ||, so the convergence of the HTE may be legitimately questioned. 
An important point to take into account in this respect is that this is true mainly for a 
light Higgs, where the approximation of keeping M w = and M H ^ manifestly fails. 
The result of [p4| |, which was based on such approximation, becomes meaningless in this 



regime, and no hierarchy among the first and the second term of the asymptotic expansion 
should be expected. Moreover, we know from Ref. j| that the size of the 0(g 4 M? /M^) 
contributions to sin 2 9^Jj and M w depends very strongly on the scheme adopted. This 
tells us that reducible contributions (products of one-loop integrals) are very important 
there. Since we know that the HTE works very well at one-loop level where the leading 
quadratic term is dominant, the reducible contribution is well approximated by the first 
two terms of the HTE. Concerning the irreducible contributions, the present two-loop 
calculations are based on two-point functions only. Unlike the case of three and four 



point-functions [29], the HTE appears to work quite well for self-energies, as has been 



demonstrated in the case of QCD corrections up to three loops in the second paper of 



S 



Ref. |3q| . Ultimately, however, the most convincing evidence on the quality of the NLO 
HTE approximation comes a posteriori from the sharp reduction of the scheme and scale 
dependence of the predictions for M w , sin 2 #e//> an d, as wm be seen later, IV 

A conclusive test of our approximation can only come from a comparison of the HTE 
with a two-loop numerical calculation. A first partial test of the HTE [^TJ for the two- 
loop electroweak corrections has been obtained by comparing the results of |3],|4|] with the 
calculation of Ref. |3D| by Bauberger and Weiglein. They calculated the two-loop self- 
energies which contain the Higgs boson together with fermions through a direct numerical 
evaluation of Feynman diagrams, without any heavy mass expansion. Although this is 
only a subset of the diagrams involving the Higgs boson, the result can be used as a 
partial test of the HTE employed here, because it includes all the diagrams containing 
both top and Higgs. The Higgs mass dependences of the predictions for M w and sin 2 O^fl 



according to Refs. M and [30] have been accurately compared in H, 31]. Very small 



deviations appear to be attributable to the use of the HTE, of the order of 1 MeV and 
10~ 5 for M w and sin 2 #e//> respectively. Such deviations are well within the range of 
theoretical errors estimated from scheme and scale dependence H. H . 



4 The Z° decay partial widths 

This section contains the main result of the paper, namely the 0(g 4 M 2 /M^) corrections to 
the partial decay widths Z — > ff, Tf, of the Z° boson. The only case which is not covered 
is the one of final 6-quarks, because it involves specific 0{g 4 M 2 / M"^) vertex corrections. 
As in our previous calculations , we first derive the result in the MS framework of Ref. 
H and then discuss the scheme dependence by considering the analogous contributions 
in two implementations of the OS scheme. 

We recall that the formulation of Ref. |J employs MS couplings evaluated at the scale 
M z while using OS masses for the physical particles. In particular, the coupling constants 
of the SU(2) and U(l) groups are written in terms of e 2 = e 2 /(l + 2 5e/e)^, the MS 
electromagnetic coupling, and sm 2 9 w (M z ) = s 2 , the MS electroweak mixing parameter. 
Even though the top quark is heavier than the Z° mass scale, it is not decoupled in the 
definition of the MS couplings, according to the convention of Ref. ||. The subscript 
MS in the charge renormalization counterterm (2Se/e) indicates the pole subtraction, a 
notation that will also be used henceforth. Although the t'Hooft scale /i is conventionally 
set equal to M z , in our calculations we leave it unspecified. This allows us to study the 
scale dependence of the results. In this framework, M w and the renormalized parameter 
s 2 are calculated from G^, a, and M z through the two relations (c 2 = 1 — s 2 ) 

K = C 2 §2 = 7TQ 

Mf 1 - ||Ap' V2G„MI{1 - Arv) ' 

The 0(g 4 M 2 /M^) contributions to the radiative corrections f^Ap and Ar w were discussed 
in Refs. In appendix B we report the explicit expression for the 0(g 4 M 2 /M^) 
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contribution to (2Se/e)^. We notice that there is a technical advantage in using MS 
couplings. Indeed, the associated one-loop counterterms, 5^e and 5^s 2 , are pure poles 
that do not contain any M 2 dependence and this reduces the number of counterterm 
diagrams to be considered. 

The amplitude for the decay of a Z° boson into two fermions / can be expressed in 
the MS scheme as 

M£ = ~iL[f}f(ff\ \j$ ~ s 2 kf(M 2 z ) Jt; |0>, (6) 

where J3 1 is the SU(2) neutral current and JOj is the photonic current^. The corresponding 
expression for the partial width Z — > ff is therefore 

r / = N i J^^fl {l - 4 Jj Q, § 2 Rek f + 8 Q) s%\ 2 } + 0{m)/Ml), (7) 

where N[ is the color factor for the fermion /, i| = ±1 its isospin quantum number, 
Qf its electric charge, and the appropriate QED and QCD corrections have not been 
explicitly indicated. 

The 0(g i M 2 /M^) corrections to kf(M 2 .) have been already considered in Refs. j|,|l8|]. 

In the MS scheme they are universal (with the exception of the b final state) because they 

can be entirely expressed in terms of the 7 — Z self-energy, as a consequence of the 

discussion in Sec.0. The overall factor fjf can be expressed as fjf = 1 + (e 2 /s 2 )(Ar/j. 1 ^ + 
(2) 

Afjf ), where the one-loop contribution is given by 

--,2 



MS 



The function V ff (M 2 z ) is defined in Eq. (12d) of Ref. @ and ReA^J(M 2 ) is given in the 
appendix of the same paper. 

In analogy to the case of k f , the MS two-loop 0(^ 4 M 2 /M^) contribution Arjf is the 
same for all light fermions with the exception of the bottom. Indeed, we have seen in 
Sec.|2] that one can choose a suitable renormalization of the unphysical sector (heavy mass 
procedure) which avoids the appearance of two-loop vertices at this order of the calcula- 
tion. In this case Ar)^ is given exclusively by the two-loop wave function renormalization 
factor of the Z° boson 

A^Arf) = Re4^)|_ . (9) 

We recall that the one-loop correction Af)^ and the one-loop MS coupling counterterms 

do not contain M 2 terms. Therefore Afj^ can be directly obtained from the irreducible 
diagrams contributing to the two-loop self-energy of the Z° plus the mass counterterm 
graphs. In appendix A the sum of these two contributions to A^l(q 2 ) is presented. The 

3 In Ref. the correction fjf is indicated as pff. 
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explicit expression of Ar)^ for arbitrary p can be obtained differentiating Eqs. ( |A3| . |A1| ) 
with respect to q 2 and is reported in Eq. ( |B1| ) (for the case M H <C M t ) and Eq. (|B2|) 
(M H ^> M w ) of appendix B. Numerically, Af/ 2 ) is small in all the relevant range of M H 
and M t values. For M t = 175 GeV, p = M z , and s 2 = 0.2315, Ap( 2 ) amounts to 7x 10 -5 at 
M H = 50 GeV and decreases almost monotonically for larger M H . The two expressions for 
light and heavy Higgs match very well around M H = 96 GeV, where they reach 3xl0 -5 . 
In the range 200 GeV £ M H £ 1 TeV, Afj^ is always less than 2xl0 -5 . 

We consider now the Z° partial widths in the OS scheme. To derive the relevant 
expressions we follow the same steps and employ the same notation as in Ref. Q. The 
idea is to use the known one-loop relations between the MS and OS schemes || in order 
to translate the MS quantities e and s appearing in Eq. ([^D into OS parameters. Following 
Ref. H, we write 



e 2 /s 2 =(G,/V2)8M 2 w [l-(e 2 /s 2 )f], 
s 2 /s 2 = l + {c 2 /s 2 ){e 2 /s 2 ) Ap, 

where s 2 = 1 — c 2 is an abbreviation for sin 2 6 W = 1 — M^/Mj and 



(10) 
(11) 



/ 
Ap 



(ReA l 



Re 



A, 



,(M, 

,(m; 



A v 



c 2 A> 



■ (Q))/Ml + V w + M 2 W B W 



MS 



; (M 



_/M, 

MS ' 



(12) 
(13) 



In these expressions A ww (q 2 ) and A zz (q 2 ) are the transverse W and Z° self-energies 



34 1 with e 2 /s 2 factored outQ. Similarly, a factor e 2 /s 2 has been extracted from the 



definition of Ap and /. The term V w + M^B W in Eq. (|12]) represents vertex and box 
diagram contributions to the p decay, modulo a factor e 2 /s 2 0. The quantities / and 
Ap are actually functions of s 2 . Iterating Eqs. (pT0| ) and (O) we can express e 2 /s 2 and s 2 
completely in terms of and s 2 . We note from Eq. (|TT|) that s 2 — s 2 = c 2 (e 2 /s 2 )Ap pa 
3 c 2 x t + ... where x t = G m M 2 /87t 2 a/2 and the ellipses represent subleading contributions. 
Thus, the replacement s 2 = s 2 [1 + (c 2 /s 2 ) 3x 4 + ...] in the one-loop components of f(s 2 ) 
and Ap(s 2 ) induces additional contributions of 0(g 4 M 2 /M^). 

We recall that in Ref. 0] two ways of writing the OS corrections Ar 
discussed. The two formulations are equivalent to 0(g 4 M 2 j 'M^ 



34 



were 



and k 

but include in different 
ways higher orders effects. The first one, indicated as OSI, generalizes the Consoli Hollik 
Jegerlehner formula |35|] for Ar to 0(g 4 M 2 /M^) by writing 



1 - Ar 



1 + 



25e 



MS 




S 2 y/2 



(14) 



where f(s 2 ) and Ap(s 2 ) are the OS counterparts of Eqs. (|T2"D and (|T^) and are defined in 
appendix B. Eq. (|TJ|) shows a structure for Ar that is very close to that obtained in the 

4 We have excluded a factor (A 7Z (Mj)) 2 from the Z° mass counterterm because it does not contribute 
at the order of the calculation. In the numerics, however, the square of the imaginary part is not negligible 
because of the additive contribution of all the light fermions, is gauge independent, and is enhanced by 
a factor n 2 . As in the case of the effective sine [|33t| , we have included it in our numerical studies. 
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MS framework ||. As the MS formulation automatically resums reducible contribution 
and therefore includes part of higher order corrections, Eq. ( |I^] ) provides an analogous 
resummation formula in the OS scheme. Moreover, similarly to MS, Eq. (|l4j) contains 
some residual /^-dependence at 0(g A ). 

A similar analysis can be performed for the Z° partial widths. Referring to Eq. (|7p, 
we note that the OS analogue kf of the form factor kf depends on the fermion /, unlike 
its MS counterpart. In Ref. kf has been presented only for a leptonic final state. In 
appendix B we generalize this result to the case of arbitrary fermion but the b quark. We 
still have to discuss the translation of the prefactor into the OS language. Indeed, we can 
rewrite ^§^|t)/| in Eq. (|7|) in terms of OS quantities keeping terms up to 0(g 4 M 2 /M^) 

and also try to resum some higher order effect, in analogy to the MS scheme. Using 
Eqs. fllOl ) and (jTTD , we obtain 





l + ^&Afj f (s 2 ) 




1 V2 


[Ap(s 2 ) - /(*»)) 



(15) 



where 
with 



Arj f (s 2 ) = Afjf\s z ) + Afjf>(s 2 ) (16) 

Af)?\s 2 ) = Ai 1 f{s 2 ) + Afjf^s 2 ) . (17) 

In the spirit of the OSI scheme, Eq. ( |i5|) follows closely the MS formulation, implies 
some resummation, and is actually /i-dependent at 0(g 4 ). Therefore, we adopt Eq. flT5] ) 
whenever we employ OSI to calculate s 2 . Explicitly, Afj ( j l \s 2 ) is the one-loop OS result 
with the overall coupling written in terms of G M , Afj^ 2 \s 2 ) is given by Eqs. ( plf ) and 

( P2|) in which s 2 has been replaced by s 2 , while Af]^l dd is the term induced by the shifts 
s 2 — > s 2 in the one-loop result and it is provided in appendix B. 

In the second implementation of the OS scheme introduced in Ref. Q, called OSII, 
Ar contains only the two-loop contributions proportional to M t 4 and M 2 , without any 
MS-like resummation, and it is strictly \x independent. In this case 

A, = A,-™ + A,- + N cXt ( 2 £!ap« - iV.^f ) , (18) 

where Ar^ is the original one-loop OS result of Refs. |||34|], expressed in terms of a and 
s 2 , Ar^ = (e 2 /s 2 )f^ — (c 2 / s 2 )(e 2 / s 2 )Ap( 2 \ The last term in Eq. fll8|) represents higher 
order reducible contributions induced by resummation of one-loop corrections, while Ar^ 
contains the corresponding irreducible components. In the same spirit of what has been 
done for Ar in the OSII formulation, we now consider a different way to write the prefactor 

"2 

1 f)/ 1 in the OS scheme. We expand the r.h.s of Eq. ( |T5j ) and retain only two-loop effects 
enhanced by powers of M t , obtaining a //-independent result: 



e 2 ,osn ZG^M 2 



s 2 c 2lVfl V2 



V2 



Af) f -f + Ap 
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T £ (MeV) 


M H 


OSI 


OSII 


MS 


65 


84.030 


84.028 


84.029 


100 


84.014 


84.012 


84.013 


300 


83.926 


83.922 


83.925 


600 


83.849 


83.845 


83.849 


1000 


83.791 


83.786 


83.791 



Table 1: Predicted values of 1^ in OSI, OSII, and MS. QCD corrections based on fi t - 
parametrization optimized. M t = 175GeV, a s {M z ) = 0.118, Aa had = 0.0280,M Z = 91.1863. 





r £ (MeV) 


M H 


OSI 


OSII 


MS 


65 


84.051 


84.062 


84.044 


100 


84.035 


84.045 


84.028 


300 


83.942 


83.949 


83.937 


600 


83.862 


83.867 


83.858 


1000 


83.801 


83.804 


83.797 



Table 2: As in Table [Lj but excluding 0(g 4 M 2 /Ml) corrections. 

+ ^M^N c x t (2Ap« + A^ 1} - 2/W) - N 2 x 2 t 



(19) 



The two OS approaches that we have introduced are equivalent at 0(g 4 M 2 /M^) but 
differ by 0(g 4 ) and higher order terms. Together with the MS scheme, they provide us 
with a broad spectrum of possibilities to study the scheme dependence of our results. 
Summarizing our discussion of the OS scheme, we write Eq. (|7|) as 

r / = N l ^ P I 1 - ^IfQfS 2 Rek f + 8Q 2 s 4 \k f \ 2 } + 0(m 2 /M 2 z ), (20) 

where we evaluate the various terms in two different ways: (i) in the OSI approach the 
prefactor P is given by r.h.s of Eq. (JE)|), while the OS sine is computed from a, and 
M z using Ar from Eq. ([14]) . The kf employed is accordingly given by Eq.(16) of Ref. @]. 
(ii) in the OSII scheme, instead, the result is obtained using Eq. (plj) for P, Ar given by 
Eq. (pD, and the factor k f given by Eq.(17) of Ref. 0]. 

Tables [1] and |2| present numerical results for the leptonic width of the Z°, Te, after and 
before the inclusion of the 0(g A M 2 / M^) contributions. These tables are the analogue 
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for Te of the corresponding tables of Q for M w and sin 2 9 l ^\. We have used exactly the 
same inputs of that paper and employed the same procedure in deriving the numerical 
results. We have neglected all final state mass effects. The two- loop |3(||37| and leading 



three-loop |38| QCD corrections are calculated according to the pit parametrization |4[], 
a procedure of implementing the QCD corrections in which the pole top-quark mass M t 
is expressed in terms of rht{nt) = fit, the MS-parameter, leading to sharply reduced 



QCD effects, and \it/M t is evaluated by optimization methods |j39| . We have verified 
that the conventional implementation of the QCD corrections in terms of the top pole 
mass leads to very close results, within 4 KeV at most. From Tables [I] and |2] we see 
that the scheme dependence of the results is strongly decreased by the inclusion of the 
0(g 4 Mt /M^) contributions, a pattern that was already observed for M w and sin 2 9 l *Ff in 
Ref. Q. In particular, we find a maximal variation of 5 KeV in Table [l], compared to 18 
KeV in Table [|. A reasonable estimate of the theoretical uncertainty due to uncalculated 
electroweak effects is therefore 0(5) KeV. Adding the QCD uncertainty, estimated to 
be ~ 7 KeV in ||, there is good agreement with the results in Table 10 of the second 
paper in ||0]. We also observe good agreement between this last reference and the error 
estimates for M w and sin 2 9^1 presented in In the case of the leptonic widths, 

however, the impact of these corrections (at most 34 KeV) remains in all schemes well 
below the present experimental error of 90 KeV (assuming lepton universality). We have 
also observed a reduction of the scale dependence of the Ti prediction in the MS and in 
the OSI scheme when the 0(g 4 M 2 /M^) contributions are included. 



5 Discussion 



The results of the 0(g 4 M 2 /M^) calculations, together with all relevant 0(aa s ) and 



the leading 0{aa 2 s ) QCD corrections [SS|, can be summarized by very simple interpolating 



formulas. Indeed, the relative uncertainty of most of the input parameters is quite small 
and their effect on the observables can be easily linearized, with the only exception of the 
Higgs mass. In Ref. formulae for sin 2 9^ and M w as functions of M H , M t , a s and the 
hadronic contribution to the running of the electromagnetic coupling have been presented 
for the three different electroweak schemes introduced in the previous sections, namely 
MS, OSI, and OSII. Here we apply the same procedure to the case of Ti and slightly 
update the results for M w and sin 2 9^Jj. The only differences with respect to Ref. [0] 
are that here we employ (i) the new value Gu = 1.16637 x 10~ 5 GeV -2 obtained after 



incorporation of the complete two-loop QED corrections to the muon decay ||T] and (ii) 
the new central value M z = 91.1867 GeV. These two changes in the inputs amount to a 
very small shift of the overall normalization, common to the three schemes. The analytic 
formulae are of the form 

sin 2 9 l e f f = (sin 2 9 l e f f ) + c 1 A 1 + c 2 A 2 -c 3 A 3 + c 4 A 4 , (21) 

M W = M° W - dt A l -d b A\- d 2 A 2 + d 3 A 3 - d 4 At, (22) 
T e = T° - 9l A l -g 5 Al- g 2 A 2 + g 3 A 3 - g A A 4 , (23) 
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Scheme 


(sm 2 ^) 


10% 


10 3 c 2 


10 3 c 3 


10% 


MS 


0.231513 


5.23 


9.86 


2.78 


4.5 


OSI 


0.231527 


5.19 


9.86 


2.77 


4.5 


OSII 


0.231543 


5.26 


9.86 


2.68 


4.4 



eff' 



Table 3: Values of (sin 2 

incorporate 0(g 4 M 2 /M^) corrections in the /x t -parametrization of QCD corrections [2]. 



o and q in Eq. (|2l|) for three electroweak schemes that 



Scheme 


M° 

w 


10 2 di 


10 d 2 


10 d 3 


10 2 d 4 


10 3 d 5 


MS 


80.3829 


5.79 


5.17 


5.43 


8.5 


8.0 


OSI 


80.3809 


5.73 


5.18 


5.41 


8.5 


8.0 


OSII 


80.3807 


5.81 


5.18 


5.37 


8.5 


7.8 



Table 4: Values of M° and d* in Eq. ([22]), in GeV, as in Table 



where 

100 GeV 0.0280 

A - ( Mt V 1 A - as{Mz) 1 

(Aa)ft is the five-flavor hadronic contribution to the QED vacuum-polarization function 
at q 2 = Mf , and (sin 2 l ^) , M^, and T° are (to excellent approximation) the theoretical 
results at the reference point M H = 100 GeV, (Aa) h = 0.0280, M* = 175 GeV, and 
a s (M z ) = 0.118. The values of (sin 2 9 l e J}) , T°, q (i =1-4) and d,-, ^ (j =1-5) for the 
three electroweak schemes of Ref. are given in Tables |3|-[5[ We present the coefficients 
in the case of the /i t -parametrization mentioned in the previous section, which was shown 
in Ref. 0] to give results very close to the more conventional parametrization in terms 
of the top pole mass. These formulae are very accurate for 75 M H 350 GeV with 
the other parameters in the ranges 170 <, M t <, 181 GeV, 0.0273 <, Aa h ^ 0.0287, 
0.113 a s (M z ) <> 0.123. In this case they reproduce the exact results of the calculation 
with maximal errors of ^s 2 // ~ 1 x 10~ 5 , 5M W <^ 1 MeV and ST^ 3 KeV, which are 
all very much below the experimental accuracy. Outside the above range, the deviations 
increase but remain very small for larger Higgs mass, reaching about 3 x 10~ 5 , 3 MeV, 
and 4 KeV at M H = 600 GeV for sin 2 ^, M w , F e , respectively. 

By comparing the normalized coefficients of lnM H in Eqs. ( pT|) and (22), we see that 
sin 2 O^ff is ~ 3 times more sensitive to In M H than M w , ~ 6.6 times more sensitive to Aah, 
almost 2 times more sensitive to M t and a s . Concerning Eg, table |5| shows a sensitivity to 
lnM H comparable to that of M w . Its present experimental determination, however, is less 
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Scheme 


1 1 


10 2 ^! 


10<? 2 


10<?3 


10<? 4 


10 2 <fc 


MS 


84.0112 


5.37 


4.76 


8.01 


1.13 


2.14 


OSI 


84.0122 


5.39 


4.76 


8.02 


1.12 


2.16 


OSII 


84.0099 


5.49 


4.75 


7.99 


1.12 


2.16 



Table 5: Values of T° and gi in Eq. (|2"3"|), in MeV, as in Table 0. 



precise and it is not expected to improve significantly in the foreseeable future. Despite 
the recent progresses in the measurement of M w , it is clear that most of the present 
sensitivity to M H still comes from the effective sine. However, with an improvement in 
the measurements of both M w and M t the situation may change in the future |7]-|9| • 

We now use Eqs. (|2lH23|) together with the latest data to determine M H by a x 2 
analysis. The experimental inputs we use are [[E|: sin 2 9 l e J} = 0.23157 ± 0.00018, M w = 
80.394±0.042 GeV, T e = 83.90±0.10 MeV, M t = 174.3±5.1GeV, a s (M z ) = 0.119±0.003. 
Concerning (Aa)h we adopt the value 0.02804 ± 0.00065 [43] but will comment later on 
the impact of some recent estimates [44]. The dependence on the scheme of our result 
is very mild. The MS scheme leads to the highest M H value. The ln(M H /100) is to 
good approximation normally distributed and we obtain in MS scheme lnM w /100 = 
—0.04 ± 0.64, corresponding to a central value M H = 96 GeV and to a 95% probability 
upper limit of about 270 GeV. Repeating the same analysis in the other two schemes, 
we obtain 95% upper limits that differ from the one in MS by at most —15 GeV. One 
should also consider the uncertainty due to uncalculated QCD corrections. The dominant 
QCD effect is linked to the leading 0{G ll M^) contribution to Ap, and in particular to the 
top quark mass definition. A simple way to take the related uncertainty into account is 
through a systematic shift of M t . Indeed, the conservative values for the QCD uncertainty 
given in [^, §] correspond to a variation of SM t = ±0.9 GeV. An increase of 0.9 GeV in 
the central value of the top mass moves the M H 95% upper bound by about +15 GeV. 
Thus M H <, 285 GeV could be taken as a 95% probability upper limit obtained using 
only indirect information and including the theory uncertainty. Finally, we note that a 
corresponding analysis carried out without implementing the 0((7 4 M 2 /M 2 / ) corrections 
gives a central value and upper bound for M H significantly larger, by about 40 and 90 
GeV, respectively. 

Following Ref. flrOfl, we can include in the M H analysis the information coming from 



the direct search experiments presently in operation at LEP. The combination of the 
search data, for center of mass energy up to y/s = 183 GeV, with the above result 
from precision measurements moves the central value for M H from 96 to 145 GeV, while 
the 95% probability upper limit including the theory uncertainty reaches 345 GeV. The 
uncertainty on the upper bound due to the simplified likelihood procedure adopted has 



been estimated in |10[ to be 0(5 GeV). 



Because of the strong correlation between M t and M H , apparent in Eqs. (21 -23), in a 
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global fit also observables insensitive to the Higgs boson can have an indirect effect in the 
M H determination. In particular, the prediction for Rb is almost independent of M H , but 
the present experimental value Rb = 0.21680 ± 0.00073 points towards a much lighter top 
(Mt ~ 140 ± 24 GeV) than most other data. A smaller M t central value favors a lighter 
M H and consequently leads to a reduced 95% upper bound. The inclusion of Rb in the 
fit, for example, would lower the M H central value and upper limit by approximately 10 
and 30 GeV, respectively. Our analysis, in any case, is based only on M w , sin 2 and 
T^, the three best measured observables which are also known to 0(g 4 M 2 /M^,). 

Several new estimates of the hadronic contribution to the running of a have appeared 
over the last couple of years ||44|| . These new investigations differ from the previous most 
phenomelogical analyses fG| mainly by the use of perturbative QCD (pQCD) down to 
energies of the order of a few GeV and by the treatment of old experimental data in regions 
where pQCD is not applicable. The combination of these two factors gives a result for 
(Aa)h with a drastically reduced uncertainty but, at the same time, with a lower central 
value. Since these new estimates remove a significant part of the theoretical error in the 
prediction of sin 2 O^ff, they further enhance its role in the indirect M H determination. As 
can be seen from the sign of the coefficient C2 in Eq. (f2~Tl), a smaller central value for (Aa)h 
prefers a heavier Higgs. There is therefore a compensation between the reduced error and 
the lower central value so that only the very high confidence level values are significantly 
reduced. As an illustration, we perform the M H analysis using the Aah determination 
with the smallest error, due to Davier and Hocker, i.e. Aa^ = 0.02770 ± 0.00016. The 
MS fit, without including theoretical errors and the information from direct searches, has 
a central value at M H = 122 GeV and an upper bound only 10 GeV smaller than the one 
obtained with the Aa^ determination of [[43|] . 



6 Summary 

The 0(g 4 M 2 / 'M^) corrections have been calculated for many accurately measured pre- 



cision observables and are now implemented in several global fit programs [45j-[47]]. The 
only important quantities whose predictions are not known at the same order are those 
connected with b- final states because their 0(g 4 M 2 /M^) contributions involve specific 
two-loop vertex corrections. 

In this paper we have presented the results of the calculation of the partial widths Ff 
(/ 7^ b) of the Z° boson to 0(g 4 M 2 /M^) in three different electroweak renormalization 
schemes. In the case of 1^, the size of these contributions is modest (up to 34 KeV) 
when compared to the present experimental error (90 KeV), but, like for M w and the 
effective sine, their inclusion sharply reduces the scheme dependence of the theoretical 
predictions. The situation is very similar for the partial widths into quarks. Concerning 
T u and T d , in the case of a light Higgs of 100 GeV the 0(g 4 M 2 /M^) corrections induce 
the variations Ar u = —63 KeV and AI^ = —85 KeV in the MS scheme, Ar u = —62 
KeV and AI^ = —73 KeV in OSI, while the shifts in OSII are approximately twice 
larger, Ar u ~ AI^ = —175 KeV. It is interesting to compare these shifts with the ones 
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due to non-factorizable 0(aa s ) corrections to the hadronic width: AT U = — 124 KeV, 
Ar^ = —173 KeV [pD| . Adding the contributions of the u,d,s,c quarks, the overall 
shift to the total hadronic width I\ may reach 0.7 MeV, to be compared to the present 
experimental error of 2.3 MeV. 

We have also described in detail several aspects of the 0(g 4 M 2 /M^) calculations, 
which had not been reported in previous publications. Furthermore, in appendix A we 
have presented explicit expressions for the two-loop 0(g A Mf /M^) contribution to the 
self-energies of the W and Z° bosons. 

The analytic results of the calculations for M w , sin 2 6*^, and 1^, together with the 
relevant QCD corrections, have been incorporated in a numerical code, from which we have 
derived simple interpolating formulae following Ref. 0. They reproduce very accurately 
the results of the code for a large range of values of the input parameters. At present, 
one of the most interesting application of high precision investigations is the attempt to 
constrain the Higgs boson mass from the data. As has been repeatedly observed 0,0-0], 
our 0(g 4 M 2 /M^) calculations do have a significant impact on the Higgs mass bound, 
because of an enhancement of the screening operated by the higher order contributions 
on the leading one-loop correction. Using the interpolating formulae, it is very simple 
to obtain an indirect determination of M H based on the latest preliminary data. The 
result of our analysis using only the indirect information from precision measurements 
gives M H < 285 at 95% C.L., where the theoretical uncertainty is also taken into account. 
The inclusion of the direct search information up to energies ^/s = 183 GeV according 
to the procedure outlined in |TD[ increases the 95% upper limit by ~ 60 GeV. As a 



final remark, we notice that the Higgs determination is very sensitive to the value of the 
effective sine employed in the analysis and that the above results have been obtained 
using the SLD+LEP average. Although recently the discrepancy between the LEP and 
SLD sin 2 9 l ^Jj determinations has somewhat receded, it should be kept in mind that their 
agreement is still far from satisfactory. 
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Appendix A 

In this appendix we present some of the intermediate results. The analytic expressions for 
the irreducible two-point functions of the W and the Z° bosons through 0(g 4 M 2 /M^) 
are extremely long and cannot be reduced to a compact form. However, after on-shell 
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renormalization of the masses inside the loops, i.e. after adding the mass counterterm 
diagrams, they simplify considerably. 

The expressions for the two-loop self-energies that follow are calculated in the 't Hooft- 
Feynman gauge and include all the irreducible unrenormalized diagrams that contribute to 
the order of the calculation plus the corresponding contributions from the on-shell mass 
counterterms of the particles in the loops. Neither gauge coupling nor wave-function 
renormalization are carried out. The renormalization of the unphysical Z° and W sectors 
is performed according to the heavy mass procedure described in Sec. ^|, so that all two- 
loop vertex and box contributions proportional to powers of M t are canceled, and only 
two- loop self-energies have to be taken into account. 

As the couplings of the SM sometimes contain masses, which we renormalize in the 
same way as the ones that appear in the propagators, a possible ambiguity arises in the 
MS scheme, because it is always possible to rewrite M w as M z c in the couplings of a 
one- loop diagram; the counterterm associated to the two forms differ in the finite part. In 
order to avoid any possible source of confusion, it is important to stress here that all our 
two-loop expressions are therefore consistent only with the choices adopted in the one-loop 
MS expressions reported in |]5,32]. Different choices would induce different 0(g 4 M? /M^) 



terms. 



We first introduce some functions which appear in the following expressions of A ww (q 2 ) 
and A zz (q 2 ) and in the formulae of the next appendix. We define: 



9\ x ) 



a/4- x (ir — 2 arcsin Jx/A) < x < 4 



2 x fx~/4^l In ( i-VEgE ) x > 4 



A{x) 



.^.g( X ) + ^yjA/x-l 0<X<4 



where we have indicated the dilogarithmic function as Li 2 (x) 
introduce 



J*dt 



ln(l-t) 



and we 



4 VIZi Ch(2 arcsin y/z) 



4Li 2 (i^) + 2 ln 2 (i^) - ln 2 (4^) + tt 2 /3 



< z < 1 
z > 1 , 



where Cl2(x) = ImL^e") is the Clausen function and A = Jl — ~. 

The function B [q 2 , m 2 , m 2 .} is defined through the one-loop integral (e = (4 — n)/2) 

d n k 1 1 



i M 2t e^ e 



+ B Q [q 2 ,m{ } mi] + 0(e) 



(Al) 
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whose analytic form is well known (see for example Ref. 
that the 0(e) part of one- loop integrals like the one in Eqs. 
two-loop expressions of the self-energies given below. 

We also introduce the following short-hand notation: 



It is interesting to note 
If) cancel exactly in the 



w 



Ml; 



w t = f 



t ■ 



h = Ml- t = M 2 ; 



t > 



As no coupling renormalization has been carried out in the following self-energies, we 
express the results in terms of unspecified sine and cosine s, c of the Weinberg angle 
which should not be understood as on-shell quantities. 

• A zz (q 2 ) 

In order to fix completely the one-loop contribution we report here the coefficient of the 
1/e pole term that is not explicitly written in Ref. || 



a (l) pole _ 
ZZ ' (167T 2 ) 



41 41 



6 c 2 



H -M f 2 + 4 

3 / 2 c 2 ' 



Ml # 



(A2) 



where M t , M w and M z are understood as OS masses. Concerning the HTE of the two- 
loop self-energy, as explained in Sec. || we have to distinguish between the case in which 
the Higgs is light relative to the top quark, and the one in which it is comparable to it. 
In the first case, the result, expressed in units N c g*/c 2 M 2 /(16tt 2 ) 2 (yu/M t ) 4e (JV C = 3), 
reads 
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Here and henceforth S t is the 0(e) part of the top mass counterterm calculated at = M t 
which cancels out in the physical amplitudes we have considered. In the case of a heavy 
Higgs boson, the HTE takes the form, in the same units: 
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The coefficient of one-loop pole is 

,2 



^(l)poie 



.9 



(16tt 2 



5 
6' 



--g 2 + ^M 2 -^M^ + Mjc 2 



(A5) 



In the case of a light Higgs, the HTE for the two-loop W self-energy takes the form, in 
units ^ 4 M 2 /(16tt 2 ) 2 (n/M t ) 4e , 
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while in the case of a heavy Higgs, we find in the same units: 
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The other self-energies used in the various 0(g A M 2 / M"^) calculations are the photon two- 
point function at g 2 = and the 7 — Z° mixing amplitude at q 2 = M 2 . They are the 
only contributions to <^e/ e lMis an< ^ ^k^ 2 \ respectively, and therefore can be read directly 
from Eqs. (|B|) and flE^). 

The limit for g 2 — > of the previous expressions can be easily obtained using the 
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expansion of B [q , a, b] in powers of g : 



alna-blnb 2 (a 2 - & 2 - 2 a Mn f ) 



S [g 2 , a, 6] = 1 + q 2 ± -j— ^ + 0(g 4 ) (A8) 

a — 2 (a — o) 

As an example of the use of the above self-energies, we now calculate the MS coun- 
terfeit! 5s 2 at 0{g 4 M 2 jM 2 w ). We recall that, following Eq. (12b) of Ref. §, 5s 2 = 
—M^/M 2 Y\ po i e , where Y is the expression appearing on the r.h.s of Eq. (|13D , before 
the poles are subtracted. Inserting the poles of Eqs. (|A3| ) and ( |A~6| ) or Eqs. ( |A"4T ) and ( |A~7| ) 
into the previous relation is not sufficient to obtain the correct result in the pure MS 
scheme because (i) no one-loop coupling renormalization has yet been carried out and 
(ii) the two-loop self-energies have been obtained using OS masses instead of a pure MS 
subtraction. The MS renormalization of gauge couplings does not introduce M 2 terms 
and at 0(g 4 M 2 /M^) it is only needed for the leading quadratic contribution. Working 
as usual in n dimensions, we find 

_ W y WM» = _ 1» _M_ML (l + l_ ln ML + 0(e) ) (A9) 
j2 r MS 6 (16x 2 ) 2 4Af2. \e 2 /I 2 + K ) j ' K > 

In order to re-express the poles of the two-loop self-energies in terms of MS quantities we 
take Eqs. fl5|) and (|Aj), insert them in M^/M 2 Y, replace M^/M 2 by c 2 (l + fyApj^) 
in n dimensions, and expand at 0(g 4 M 2 / 'M^) ■ The extra term at this order reads 

(2) N c g 4 M 2 (I 1 , M 2 n . \ ( „ 13 1 \ , . 



Combining the poles of the two-loop self-energies in Y together with Eqs. (|A9|) and Y t 
we find the correct MS counterterm 
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which can be checked with ||49|| . Combining the finite parts of these three components, 
instead, one obtains Ap^ given in ||. 



Appendix B 

In this Appendix we present the two-loop 0(g 4 M 2 /M^) contributions to the electroweak 
corrections entering the Z decay partial widths. 

We begin by considering the f) and k corrections of the MS framework. We have for 
M H < M t in units iV c /(167r 2 ) {a/Ans 2 ) M 2 /{M 2 Z &) 
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In the case of M H ^> M z we obtain in the same units 
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The correction can be written as kf = 1 + (e 2 /s 2 )(Akf 



A (21 

Akj ) where the two- 



loop part has been given already in Ref. M. We report here again for completeness. 
Considering M H > M t in units A^ c /(167r 2 )(d/(47rs 2 ))M 2 /(c 2 M|) one has 
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In the case of a light Higgs mass, this becomes 
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Numerically, Re AM 2 )(M|) is tiny, and the use of the expansion for light Higgs does not 
lead to any appreciable difference over the other one. 
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We give also the two-loop finite part of the electric charge counterterm, which is needed 
to calculate the MS electromagnetic coupling. We have in units N c a 2 / (16tc 2 s 2 ) M 2 / (Mjc 2 ) 

25 { - 2 \ 
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that in the case of a light Higgs reduces to || — 16 — j| In ^f-. 

We proceed to the discussion of the corrections of the OS scheme. Recalling Eq. ([1 



the additional two-loop contribution in Arj^ is given by (units iV c £t/(167r 2 ): 
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where V add contains the shifts of the one-loop vertices. 
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and the functions fv(x) and gv(x) are defined in Eqs.(6d) and (6e) of p2|. 

The form factor kf was considered in Ref. @] only for the leptonic case. The general- 
ization of Eq. (A4) of that reference to arbitrary fermions is (units iV c £ t /(167r 2 )) 

M% =^§^ + 8c 4 - 2 c 2 VA^l (3 + 4 c 2 ) arctan(- ? =A= T ) - fc 2 In z t 

l_ 6/ 3g +8 g2_ 8c 4g2 17 2 
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+c 2 (l-Q / / / 3 )Mc- 2 )-ynr. (B8) 

To complete the discussion in the OS framework we need also the corrections Ap(s 2 ) 
and f(s 2 ). They have been considered in Ref. @]. Recalling the definitions 

(B9) 

Ap(s 2 ) = Ap«(s 2 ) + Ap( 2 )(s 2 ), (BIO) 

where f^(s 2 ) = f^(s 2 ) + ffj d (s 2 ), Ap^(s 2 ) = Ap( 2 \s 2 ) + Ap a 2 J d (s 2 ), we have that 
(<V V2) 8Af£ Ap( 2 ) is given by Eqs.(10a,b) of Ref. § multiplied by N c x 2 while Apfi(s 2 ) 
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is reported in Eq.(A2) of Ref. HQ. The amplitude (e 2 /s 2 )f^ is given by the sum of 
(25^e/e) given in Eq. QB|) and Eqs.(7a,b)Q of Ref. § multiplied by (a/7rs 2 )V c x t . The 
term /i^(s 2 ) is given in Eq.(Al) of Ref. 0. 
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